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Abstract. Recent far-infrared mapping of mass-losing stars by the AKARI 
Infrared Astronomy Satellite and Spitzer Space Telescope have suggested that 
far-infrared bow shock structures are probably ubiquitous around these mass- 
losing stars, especially when these stars have high proper motion. Higher spatial 
resolution data of such far-infrared bow shocks now allow detailed fitting to 
yield the orientation of the bow shock cone with respect to the heliocentric 
space motion vector of the central star, using the analytical solution for these 
bow shocks under the assumption of momentum conservation across a physically 
thin interface between the stellar winds and interstellar medium (ISM). This 
fitting analysis of the observed bow shock structure would enable determination 
of the ambient ISM flow vector, founding a new technique to probe the 3-D 
ISM dynamics that are local to these interacting systems. In this review, we 
will demonstrate this new technique for three particular cases, Betelgeuse, R 
Hydrae, and R Cassiopeiae. 



1. Interactions between Stellar Winds and the ISM 

Most stars lose their surface matte r into the surrounding interstellar space in one 
way or another via mass loss (e.g. IWillsorJl2000l : iKudritzki & PuldteOOOT ). Such 



mass loss processes intimately relate these massdosing stars to the interstellar 
medium (ISM), which these stars enrich with ashes of nuclear burning taking 
place inside these stars. Indeed, mass loss is the critical proces s that one must 



consider when it com es to the mass budget of the ISM (e.g. ISedlmavrl 11994 ; 
iMatsuura et alJl2009ft . 



When we think of such ISM enrichment by stellar mass loss ejecta, however, 
we often naively expect that stellar ejecta will eventually meet with the ISM and 
get merged into it without giving serious thought about the process of the stellar 
ejecta merging with the ISM. On the first order approximation, the wind-ISM 
merger occurs where the ram pressure of the stellar wind is balanced by the 
pressure of the ambient ISM: a shock is expected to form if the wind velocity 
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is greater than the ambient sound speed (t> w > v c ) and a pile-up of the ejecta 
is anticipated otherwise (t> w < v c ). Either way, a density enhancement ensues 
at the interface between the stellar wind and ISM. Such density enhancements 
can be detected if sufficient amount of radiation is induced within the wind-ISM 
interface regions. Since the wind-ISM interface often occurs surrounding the 
mass- losing star, we typically detect such density enhancements in the form of 
concentric arcs around the central star. 

Such arc-shaped circumstellar structures had already been found around lu- 

minous OB stars and Wolf-Rayet a nd red supergiant stars (jStencel. Pesce. h Hagen Bauer 
19881 : IVan Buren k McCravlll988l ) in the far-infrared (far-IR) All-Sky maps ob- 



tained by the Infrared Astronomy Satellite (IRAS). Far-IR emission from these 
arcs is probably mainly due to thermal emission of the cold dust component of 
the density enhancement at the wind-ISM interface whose shock-heated temper- 
ature peaks at far-IR. There may be contributions from low-excitation atomic 
lines such as [O I] 63/iin, [O I] 145;um and [C II] 158/im at these wavelengths. 
However, the exact emission mechanism of these far-IR bow shocks remains 
unclear. While an attempt to identify their spectr oscopic nature is c urrently 
on-going with the Spitzer Space Telescope (Spitzer ; IWerner et al.l 120041 ) . this is 
beyond the scope of the present paper. Below, we will restrict our discussion to 
concentrate on the structure of these stellar wind-ISM bow shocks and what we 
can learn from it. 



2. 3-D Motions of the Star, Shock, and the ISM 



IRAS discoveries of these bow shock structures at the wind-ISM interface were 
certainly astounding. Nevertheless, the spatial resolution of the IRAS observa- 
tions (typically 2' to 5') was not sufficient enough to allow detailed structural 
analyses of the far-IR shock surfaces. To pursue this avenue of research, we had 
to wait until the far-IR renaissance of this decade wit h the coming of Spitze r 
and the AKARI Infrared Astronomy Satellite ( AKARI: iMurakami et aill2007l ). 
which permit sub-arcminute spatial resolution in the far-IR. 

This far-IR renaissance in context of the present study began with a Spitzer 
discovery of a stellar wind bow shock arc around R Hydrae (R Hya), a Mira-type 
evolved star ( Ueta et"ai1l2006l . Figure [T]), followed by a revisit to Betelgeuse, one 
of the first IRAS discoveries of a far-IR stellar wi nd bow shock, usi ng the Far- 
IR S urveyor instrument (FIS; iKawada et al.ll2007l ) aboard AKARI (jUeta et al 



2008, Figure [2]) . New far-IR images of these sources clearly reveal the character 
istic arc-shaped structure of a bow shock at the interface between stellar winds 
and the ISM at sub-arcminute spatial resolution. For these bow shocks to occur, 
the mass-losing central star has to be moving relative to the ISM that is local to 
these objects. Indeed, these stars exhibit proper motion that is consistent with 
the orientation of the bow shocks (i.e. toward the apex of the bow) 



Meanw hile, theoretical stu dies have been done both analytically (IWilkinl 

19961 . 12000T) and numerically (e.g . iMac Low et al.lll99ll;lDgani. Van Buren. fc Noriega-Crespo 
1996; Bl ondin fc Koerwerl TT998I : IWareing. Zijlstra. fc O'Brienl 120071 ) to under- 
stand the structure of these bow shocks. Remarkably, if one assumes momentum 
conservation across a physically thin (i.e. radiative-cooling dominating) shock 
layer, one can express the shape of the bow shock analytically as a function of 
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Figure 1. Spitzer maps of R Hya at 70 and 160/im (top left and middle, 
respectively; also shown are the point-response function at these wavelengths 
at the botto m right) along with an Infrare d Space observatory map at 60/^m 
(bottom left; Hashimoto & Izumiura 1997) and an IRAS map at lOO^m (bot- 
tom middle; HiRes -processed) as well as a SHASSA Ha map (top right; 
iGaustad et al.ll2~001[ ). Tick marks indicate angular offsets in arcseconds with 
respect to the position of the star. Surface brightness in MJy sr -1 is indi- 
cated by the scale at the bottom of each panel. A parabolic curve, which 
roughly represents the shape of the bow shock, is displayed by the dashed 
lines. The arrow at the bottom left corner shows the dir ection of the prope r 
motion of the sta r, (iig^ns) = (— 57 .68, 12.86) mas yr _1 ([van Leeuwenll2007T) . 
Reproduced from IUetaet~aT1 (|2006h . 



the latitudinal angle 9 measured from the dir ection of the b ow apex with respect 
to the position of the central star as follows (Wilkin 19961 ): 



R{9) = R 



V3(l — cot i 



Sill! 



(1) 



where Rq is the stand-off distance between the star and bow apex denned as 



Rri 



4vrp IS Mf* 



(2) 



for which M is the rate of mass loss, v w is an isotropic stellar wind velocity, Pism 
is the density of the ambient ISM, and v* is the space velocity of the star. 

The observed shape of the bow shoc k is really a p rojection of the bow shock 
cone (eqn{T|) onto the plane of the sky ( Wilkin 1996). For instance, while the 
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Figure 2. AKARI/FIS maps of Betelgeuse in the SW bands - N60 (65/im) 
and WIDE-S (90/im) at 15" pixel" 1 scale - and in the LW bands - WIDE-L 
(140/xm) and N160 (160/im) at 30" pixel -1 scale - from left to right, re- 
spectively. Tick marks indicate angular offsets in arcminutes with respect 
to the position of the star. Surface brightness in MJy sr _1 is indicated by 
the scale at the top of each panel. The revised Hippar cos proper motion is 
(UniUs) = (27.54 , 11.30) mas yr _1 (jvan Leeuwenll200"7h . Reproduced from 
lUeta et al.1 (|2008H . 



R Hya bow appears very close to an edge-on bow shock cone (Figure Qj) the 
Betelgeuse bow looks rather circular (Figure [2]) , suggesting an inclination angle 
far from edge-on. At far-IR, the surface brightness is proportional to the dust 
column density along the line of sight. If the bow shock is physically thin we 
can expect that the column density becomes the highest where the bow shock 
layer intersects with the plane of s ky including the centr al star, based on dust 
radiative transfer calculations (e.g. lUeta Sz Meixner 2003). Hence, one can con- 
fidently determine the orientation of the bow shock cone by fitting the apparent 
bow shock shape with the Wilkin bow shock cone equation (eqn(T]), taking into 
account the inclination angle of the bow. This Wilkin solution fitting, therefore, 
would yield the heliocentric orientation of the bow shock cone with respect to 
the observer. This orientation of the bow represents a relative motion between 
the mass-losing central star and the ISM local to the star. At the same time, the 
central star's heliocentric space motion can be defined observationally via mea- 
surements of its proper-motion and radial velocity. Naturally, the probability of 
detecting a bow shock is increased if the star has a large space motion. 
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Here, we have two vectors of motion concerning a mass-losing star and its 
ambient ISM. One is the heliocentric relative motion of the star with respect to 
the ambient ISM and the other is the star's heliocentric space motion. These two 
vectors are not necessarily the same because the ambient ISM can move/flow 
into a certain direction irrespective of the star's space motion. Therefore, if 
one compares the result of the Wilkin solution fitting of the bow shock shape 
(i.e. the heliocentric relative motion of the mass-losing central star with the 
ambient ISM) and the observed space motion of the mass-losing central star 
(i.e. the heliocentric space motion of the star), one can deduce the heliocentric 
flow vector of the ambient ISM in principle. 



3. Wilkin Solution Fitting 



Now, we would like to review the potential of this investigation into the helio- 
centric flow vector of the ambient ISM, wh ich is local to a mass-losing moving 
star, based on the Wilkin solution fitting ( Wilkini 1 19961 ) to the apparent bow 



shock shape for the following cases of three particular stars. 
3.1. Betelgeuse (Red Supergiant) 

Using the AKARI/FIS image at the most sensitive WIDE-S (90/xm) band (sec- 
ond from left in Figure [2]) , the best-fit of the Wilkin solution fitting yields the 
inclination angle of the bow shock cone with respect to the plane of the sky 
(#inci) of ±(56° ± 4°), position angle of 55° ± 2° (east of north), and stand-off 
distance of 4f8 ± O'.l. The leading ± sign of 8[ nc \ indicates the unresolvable de- 
generacy of the fitting, i.e., the bow shock cone points either away from or to 
us but the fitting is unable to determine which. Nevertheless, the Wilkin fitting 
defines a space motion vector of the star relative to the ambient ISM (albeit the 
degeneracy), which we define as v*ism- 

Meanwhile, using the most recent, multi-epoch prope r-motion study of 
Betelgeuse using VLA data. lHarper. Brown. Sz Guinanl (|2008l ) derived (v a , v$) = 
(23.3,8.9) km s" 1 at 197pc. With the average radial velocity of the star (v ra d) 
at 20.4 ±0.4 km s _1 (pointed away from us), these numbers yield another space 
motion vector of the star having 6*i nc i = 40° and the position angle of 69°. This 
space motion vector is of course relative to the Sun (i.e., the heliocentric correc- 
tion applied), which we denote as v*©. Using these two space motion vectors of 
the star in two distinct frames, we can then derive the heliocentric flow vector 
of the ambient ISM, vism©> by 

v ISM0 = V* Q — VrfSM- (3) 



Before we proceed any further, however, we have to break the degeneracy of 
v*ism • Since Betelgeuse is in the vicinity of the Orion Nebula Complex ( ONC), 
we ha ve relatively complete information about its local environment (e.g. lO'Delll 
(2001)). One of the bases for the present study is the presumed existence of the 
ISM flow local to the mass-losing moving star. Therefore, there must exist a 
source of this ISM flow, and the star forming regions in the ONC are most likely. 
If we assume there is an overall isotropic outflow from the ONC, Betelgeuse's 
proper motion suggests that it has been traversing across this ISM flow in front 
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Figure 3. 3-D relationship among the heliocentric space motion of Betel- 
geuse (v»q), heliocentric flow of the ambient ISM (vism©), and space motion 
of Betelgeuse relative to the ambient ISM (v*ism)' The sphere at the lower 
right corner indicates the position of Betelgeuse at the origin of the equatorial 
coordinates (a, 5, p). 



of the ONC (see Figure 4 of lHarper. Brown, Guinanll2008l ) . If this is the case, 
the apex of the bow shock cone around Betelgeuse should point away from us, 
i.e., #i nc i = —(56° ± 4°). Thus, the degeneracy is lifted and relevant quantities 
are only dependent on the ambient ISM density, prSM- 

Using R n = 0.3pc (lUeta et al-I^POST ) and M = 3 x 10 -6 M yr" 1 as derived 
by lHarper. Brown. Guinanl (I2008T) at the best-est imated distance of 197pc to 
Betelgeuse and v w = 17 km s -1 (jBernat et al1ll979l ). Betelgeuse's space motion 
with respect to the ambient ISM turns out to be 



v*ism 






(km s 



(4) 



At the las t step, we adopt the ISM numbe r density in front of the ONC, raisM = 
0.3 cm -3 ( Frisch. Sembach. &: Yorklll990n . Therefore, the heliocentric ambient 
ISM flow at Betelgeuse is 



v ISM0 = V* — VrfSM 






(km s" 1 ) (5) 



with 11*0 = 32.4 km s , f*isM = 72.1 km s _1 , and fisM© = 42.7 km s _1 . 3-D 
relationship among these three vectors is shown in Figure O Hence, Betelgeuse 
is traversing across the ambient ISM flow (of 42.7 km s _1 ) at the heliocentric 
space velocity of 32.4 km s _1 , resulting in the space velocity relative to the 
ambient ISM at 72.1km s" 1 . In the frame of stellar winds, therefore, the flow 
velocity of the ambient ISM is 89.1 km s _1 , which can induce a strong shock. 
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Betelgeuse is located about 200 pc away from the ONC. At 89.1 km s _1 
the crossing time is roughly 2.2 Myr, which is consistent with the crossing time 
of Betelgeuse in front of the ONC. Thus, the assumption of an outflow from 
the ONC being responsible for the stellar wind-ISM bow shock of Betelgeuse is 
indeed plausible in retrospect. 



3.2. R Hydrae (AGB Star) 

Using the Spitzer 70fim map (top left in Figure [1]), the best-fit of the Wilkin 
solution fitting yields 9- mc \ = ±(59° ± 4°) and the position angle of —62° ± 3°, 
defining v^ism- The stand- off distance of lf6 zfc O'.l translates to O.lpc at the 
adopted distance of 165pc (jZiilstra. Bedding. Sz Mattel I2002T) . The updated 
Hipparcos proper motion measurements by Ivan Leeuwenl ( 20071 ) give (v a ,vg) = 
(—45.1,10.1) km s -1 at 165pc. With v rad = —10.4 km s _1 (pointed to us), 
another space motion vector has 6>i nc i = 13° and the position angle of —77°. 

Recent CO observations have revealed a lopsided emission profile, indicat- 
i ng that there is a l arger amount of receding matter than approaching matter 
( Tevssier et al. 20061 ) . In context of a bow shock cone, this means that the cone 
is pointed away from us (hence more of the receding component). With this the 
degeneracy of v*ism is removed and we derive 



^rad 

VJSM ■ — j Va 




(km s 



(6) 




Figure 4. [Left] The vismq-^ism plot for R Hya showing behavior of vism© 
as a function of nrsM- [Right] 3-D relationship among v*© vism© and v*ism 
for R Hya and its ambient ISM. The sphere at the lower left corner indicates 
the position of R Hya at the origin of the equatorial coordinates (a,d, p). 
Depending on the choice of the unknown tiism involved in v^ism, the solution 
to vism© differs. Displayed are (1) nisM = 20 cm -3 (thick solid line, v*ism = 
7.8 km s _1 ), (2) n ISM = 5 cm~ 3 (solid line, v rfSM = 14.2 km s _1 ), and (3) 
«ism = 0.5 cm -3 (dashed line, v„ IS m = 49.0 km s _1 ). 
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as a function of the remaining unknown, nisM- Unfortunately there is no definite 
observational diagnostic for the value of wism- 

However, the vismq-^ism plot (Figure HI left panel) shows that the mini- 
mum value of vism© is 45.0 km s -1 when tiism = 6 cm -3 and that uism© remains 
roughly ~ 45.0 km s _1 for larger tt-ism but increases for smaller nisM- 3-D re- 
lationship among these three vectors is shown for three different values of nisM 
(Figure [U right panel). The space velocity of R Hya relative to the ambient 
ISM is 34.9/^/nisM km s -1 , and hence, the smaller nisM becomes the larger 
t>*TKM becomes. In the frame of stellar winds (v w = 10 km s _1 , iKnapp et al.l 



19981 ). therefore, the flow velocity of the ambient ISM with respect to R Hya is 
10 + 34.9/^/nisM kni s _1 ) which can reach ~ 50 km s _1 for nisM < 0.8 cm -3 . 
Typically lower nisM values are preferred for high galactic latitude objects such 
as R Hya. Thus, the i>ism© value (and hence, the ISM flow velocity relative to 
R Hya) can be even higher. The next step in this line of research is to estab- 
lish shock diagnostics for low density, low velocity shocks and obtain the shock 
velocity in order to fix nisM- This will pre sent a new techn i que t o probe local 
ISM, complementing previous work such as lLallement et ail ( 20031 ). 



3.3. R Cassiopeiae (AGB Star) 

Far-IR images of R Cassiopeiae (R Cas) show an elliptically elongated shell 
around the central star that is located off-cente r of the shell, bu t the shell does 
not show any obvious signature of a bow shock ( Ueta et"aill2009l ). However, the 



wind crossing time for the shell is much longer than the time it took for the 
star to move to the current off-center position from t he center of the she ll at the 
velocity known from proper motion measurements (Ivan Leeuwedl2007l ). Thus, 



we speculate that the shell shape is maintained by the interactions between the 
stellar winds and ambient ISM and that the absence of an apparent bow shock 
is due to its inclination angle. 

The best-fit of the Wilkin solution fitting yields 9 ind = ±(68° ± 2°) and the 

position angle of 74° ± 2°, defining v*ism- The stand-off distance of l'A ± Oil 

trans lates to O.lpc at the adopted VLBI-measured distance of 176p c (IVlemmings, van Langevelde. &: Diamc 
2005). The updated Hipparcos proper motion measurements by Ivan Leeuwen 



(2007) give (v Q ,v 5 ) = (71.4, 14.6) km s" 1 at 176pc. With v Iad = -22.9 km s" 1 



(pointed to us), another space motion vector has #i nc i = 17.5° and the position 
angle of 78.4°. Keeping the degeneracy of v+ism we derive 



visM© = v a = 71.4 - -— = 7.9 (km s" 1 ) (7) 



as a function of the remaining unknown, nisM- The vismq-'iism plot in Figure 
[5] shows that t»isM© reaches the minimum of 59.5 km s" 1 at nisM = 0.22 cm~ 3 
when the bow shock cone points to us or of 76.2 km s _1 at uism = 18 cm~ 3 
when the bow shock cone points away from us. In the frame of stellar winds 
at 12 km s" 1 (|Knapp et alJ ll998 N l. the space velocity of R Cas relative to the 





ambient ISM is 12 + 22.5/- v Mism km s 1 , which can reach 60 km s 1 if the bow 
shock cone is pointed to us and nisM = 0.22 cm -3 . For low galactic latitude 
objects such as R Cas, relatively higher nisM values are expected. Thus, we may 
expect a value close to the lower limit of 60 km s _1 . 
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4. Summary 

When a bow shock structure is observed around a mass-losing star, the central 
star's proper motion is often consistent with the direction to the apex of the 
bow. To the first order, this gives an intuitive picture of the mass-losing star 
moving in the otherwise stationary ambient ISM. Higher spatial resolution far-IR 
images of such stellar wind bow shoc ks now allow d etailed fitting of the structure 
given the analytic shape derived by IWilkfn ([l996). The Wilkin solution fitting, 



performed for three particular cases, has demonstrated that the heliocentric 
space velocity vector of a mass-losing star and the orientation of the bow shock 
cone do not necessarily align with each other, suggesting the existence of a flow 
of the ambient ISM local to the star. 

Since the Wilkin solution fitting is inherently degenerate and other inde- 
pendent observations are needed to break this degeneracy (i.e. determination 
the orientation of the bow shock cone with respect to the plane of the sky) and 
to derive the density of the ambient ISM, the analysis is not self-sufficient. Nev- 
ertheless, three cases we have shown here unequivocally indicate that flows of 
the ambient ISM are rather prevalent and high proper-motion stars are often 
traversing such ambient ISM flows, generating bow shocks as a result. There- 
fore, stellar wind-ISM bow shocks are potentially excellent probes of the 3-D 
ISM dynamics local to these shock structures. 

Recent far-IR imaging surveys of mass-losing stars such as AKARFs own 
MLHES Mission Program (PI: I. Yamamura) and its extension for Spitzer {Spitzer- 
MLHES, PI: T. Ueta) have identified a number of sources that seem to possess 
bow shock structures. Using these data sets, one can additionally learn the 
stellar wind-ISM interactions besides the mass loss histories (from the undis- 
turbed parts of the extended dust shells) that these investigations are originally 
designed to do. With AKARTs All-Sky Survey data, there will be even more of 




Figure 5. The vism©-«ism plot for R Cas showing behavior of vism© as a 
function of nisM- The solid line is the case where the bow shock cone points 
to us and the dashed line is the case where the bow shock cone points away 
from us. 
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these far-IR bow shocks associated with high proper-motion stars. Therefore, it 
appears that large-scale investigations of far-IR bow shocks in context of prob- 
ing the 3-D ISM dynamics should be pursued at least for the solar-neighborhood 
and we intend to follow up this new avenue of research in near future with ex- 
isting AKARI and Spitzer data as well as data from near-future opportunities 
such as Herschel Space Observatory and Stratospheric Observatory for Infrared 
Astronomy. 
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